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We study a 2- and a 4-dimensional nearly integrable symplectic maps using a singular 
perturbation method. Resonance island structures in the 2- and 4- dimensional maps are 
successfully obtained. Those perturbative results are numerically confirmed. 

§1. Introduction 

Chaotic Hamiltonian systems with many degrees of freedom are not only of 
interest in view of the foundation of dynamical system, but also have a lot of ap- 
plications. Classical Hamiltonian mechanics has been applied to various physical 
problems, such as chemical reactions, plasma physics and solar systems. 1 ) Although 
it is important to study such systems, we cannot understand the global flow gen- 
erated by a Hamiltonian in the phase space due to the high-dimensionality. For a 
low-dimensional system, the phase space is visualized and then one overinterprets 
high-dimensional Hamiltonian flows beyond ones in a 2-dimensional phase space, but 
such imaginations are not guaranteed. 

On the other hand, the origin of the global Hamiltonian chaos is an overlap of two 
resonance islands in the phase space. 1 ) Then a fundamental problem for the study 
of Hamiltonian chaos is to analyze a resonance structure. To tackle this problem 
for systems with many degrees of freedom, we take a singular perturbation method. 
Using a systematic singular perturbation method, we can show the functional form 
of an island structure even for high-dimensional systems without information on the 
geometry of the phase space. 

In this article, we take a renormalization method which is one of singular pertur- 
bation methods. All sorts of renormalization methods remove secular or divergent 
terms from a naive perturbation series by renormalized integration constants of the 
non-perturbation solution, and each prescription of the method does not depend on 
the detail of a given system. 2 )~ 5 ) One of the reformulated renormalization method 
which we take here is easily applied to non-chaotic systems and also chaotic maps. 5 -* -7 ) 
Using this reformulated renormalization method, we tackle the standard map defined 
in 2-dimensional phase and the Froeschle map defined in 4-dimensional phase space. 
Both maps are chaotic and are known for the standard models to study Hamiltonian 
chaos. 
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§2. The Standard map 

To begin with, we study a resonance structure in the standard map 1 -* (x n , y n ) i— > 
(x n+1 ,y n+1 ) 

x n+l =x n +y n+l^ ( mod2 ^) 
y n+l =y n + £ sin ( x «) . ( mo d27T) 

Here x n and y n are the canonical variables, n(€ Z) denotes the discrete time and 
e(0 < £ <C 1) is the small parameter, and this map is symplectic (dx n+1 A dy n+1 = 
dx n A dy n ). Eliminating y variable, we have 

Lx n = esm(x n ), (mod2vr) (2-1) 

where Lx n := x n+1 — 2x n + x n ~ l . A naive perturbation solution shows where a res- 
onance island appears, and the structure can be obtained using our renormalization 
method. In our case, the naive perturbation series expansion x n = x(°) n + e^t 1 )" + 
£ 2 x {2)n + ( e 3) ig obtaine d D y solving Lx^ n = 0,Lx^ n = sm(x^ n ) , Lx^ n = 
x W n cos(x^ n ), • • • , where (mod 2ir) operation is taken into account in each equa- 
tion. The non-perturbative solution is given by x(°) n = a + nP, where a and P 
are the integration constants. Here P is the value of y^ n , due to the relation 
y(o)n+i ._ ^.(Ojn+l _ x (o)n _ p Tbgn the values of P and a are in [0, 1). The naive 
perturbation solution up to 0(e 2 ) is classified by the value of P as follows 

„ _ — sin(a + nP) 2 sin(2(a + nP)) 
x n = a + nP + e V- — — - + e l 



4sin 2 (P/2) 32sin 2 (P)sin 2 (P/2) 
+C(e 3 ),(mod2vr), (P^O.tt) 

x n = a + £ S _^l n 2 + £ 2^M^) (n 4 _ n2) + 0(e3)j (mod27r) ( p = Q) 

2 48 

- sin(a + Tin) 2 - sin(2a) 2 

x = a + 7rn + e ; h e n 

4 16 

+C(e 3 ), (mod2vr) (P = tt) (2-2) 

Here we concentrate on the case of P = tt as an example. In the series (|2-2p . the 
secular behavior (oc e 2 n 2 ) appears. The renormalization method removes the secular 
terms in a systematic way. First we define the renormalized variable a n as follows, 

9 — sin(2a) 9 

a n :=a + e 2 r^n 2 . 2-3 

lb 

This definition ()2-3|) should include all secular term up to 0(e 2 ) in (|2-2|) and the 
integration constant, a, coming from the non-perturbative problem. The way to 
define a renormalized variable is quite similar to one in the case of a differential 
equation. 5 ) The renormalization equation is the map satisfying a n perturbatively, 
and we impose here that the renormalization map which we would like to obtain is 
symplectic and autonomous. This is because the original map (|2T|) has these two 
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properties. Prom the definition of the renormalized variable (|2-3jl . we have 



La n = L[ a + e 



sin(2a) 
16 



-n 



,sin(2a) 



To obtain the map containing a n only, we substitute a = a n + 0(e) coming from the 
definition (|2-3f) into the relation La n = — e 2 sin(2a)/8. Then we have the following 
renormalization map, 

2 sin(2a n ) 



Lai' 



(2-4) 



When we introduce the new variable b n defined as b n+1 := a n+l — a n , it turns out 
that this renormalization map ()2-4(l is automatically symplectic (da n+1 A db n+1 = 
da n A db n ). It is noted here that the reduced map (|2-4j) is also obtained using 
another singular perturbation method, and the similarity of the two maps (|2Tj) and 
()2-4|) gives a self-similarity in the phase space. 8 ) 

From Eqs. I|2-2|) and (|2-3|) . the relations between the renormalized variable a n 
and the original variables (x n ,y n ) are obtained as 

x n = a n + nn - | sin(a n + irn) + 0(e 3 ), (mod27r) 

y n = ( a n - a™" 1 ) + 7T - | sin(a n + vrn) + 0(e 3 ). (mod27r) 

The unstable manifold of the renormalization map reproduces the resonance island 
located near y = tt in the original system (|2-lj) . Fig. ^ shows that the phase space 
reconstructed using our renormalization method is close to one obtained by the 
original map even with e = 0.98. Intuitively, this perturbative analysis is valid only 
for the case < e <C 1, however, Fig.l shows that the validity of our renormalization 
map seems to be wider than < e <C 1 practically. To determine the actual validity 
of this analysis, we need a rigorous mathematical estimate. At any rate, it turns out 
that our renormalization method is also useful even for chaotic maps, and that our 
analysis could not be restricted in 2-dimensional maps. 




Fig. 1. The resonance island near y — tv. For both figures the value of e is 0.98, and the number of 
the initial condition is 1. (a): using the renormalization method, (b): obtained by the original 
map l|2T^ . 
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§3. The Froeschle map 

Next, we would like to study a 4-dimensional symplectic map. It is expected 
that a similar analysis gives a functional form of a resonant island even in a high- 
dimensional map, and it will be confirmed later in this section. The map which 
we study here is the Froeschle map, 9 ) (x?,a;S,y5*,j#) i-> (x? +1 , x% +1 , y^ +1 , y% +1 ), 

y n+l = y n + e_M + eC_ gin(27r(x n + ^ (modl) 

f A fC 
V2 +l =V2+^ sin(27r^) + — sin(2vr(x? + a#)), (modi) 

x n+l = x n + y n+l^ (modl ) 

x n+l = x n + y n+l_ ( modl ) 

Here n(€ Z) denotes the discrete time, x^,x 2 and y™,y% are the canonical variables, 
Ai,A 2 , C are the parameters and e(0 < |e| <C 1) is the perturbation parameter. The 
Froeschle map 9 ) has numerically been studied for the Arnold diffusion 1 ) as one of 
canonical 4-dimensional maps in several papers. 10 ) Eliminating and y 2 variables, 
we have 

Lx] = ^ sin(27rx") + ^ sin(2vr(x? + x^)), (modi), (j = 1, 2) 

where Lx™ := x™ +1 — 2x™ + (j = 1,2). The naive perturbation solutions 

x l = x f^ n +ex^ n + 0(e 2 ), (j = 1, 2) are obtained by solving the following equations, 

Lxf )n = 0, (modi) (j = 1, 2) 

Lx (l)n = il gin(27r:r (0)n ) + { 0) n +4°) n )).(modl) (j = 1,2) 

J 2tt j 2tt 

The non-perturbative solutions are x^p = cij+nPj, (j = 1, 2). Here aj and Pj are the 
integration constants, and the values of them are in [0,1). The naive perturbation 
solutions are classified by the values of Pi and P 2 , the solutions up to 0{e) are 

T»-n i nP . JAj MMaj+nPj)) C sm(2ir(a l +nP 1 ) + 2iT(a 2 +nP 2 ) \ 
j a ^ + n ^ + %47r cos(27rP j ) - 1 + 4vr cos(2vr(P 1 + P 2 )) - 1 J 
+0(e 2 ), (modi), (P x + 0, P 2 / 0, Pi + P 2 / 0) (3-1) 

/ Aj sin(27r(a,- + nPj)) C 9 . .„ . .A 

Oj + nPj+e — , J r 3 — + — n 2 sm 2vr(ai + a 2 ) 

J J V 47 i" cos(2irPj) - 1 4vr v v 7 7 

+0(e 2 ), (modi) (Pi / 0, P 2 + 0, Pi + P 2 = 0) (3-2) 
aj + e^^-n 2 sm(2ir a j) + ^" n2 sin(27r(ai + a 2 ))^ 

+0(£ 2 ), (modi) (Pi = 0, P 2 = 0) (3-3) 
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and the solutions at Pi = 0, P 2 7^ are calculated as 

n f M. 2 • /« \ C sin(27rai + 27r(a 2 + nP 2 )) 

= ai + e — !-n 2 sm(2vrai) + 

\ 47T 

+0(e 2 ),(mo 
X2 = a 2 + nP 2 + e 



4"7T 47T COs(27rP 2 ) ~ 1 

+C(e 2 ),(modl) (3-4) 
A 2 sin(2vr(a 2 + nP 2 )) 



47r cos(27rP 2 ) - 1 
C si„(2 TOl+ 2,(a 2+ nP 2 )) X + 
47r cos(27rP2j — 1 / 

The solutions at Pi ^ 0, P 2 = are obtained by exchanging the suffixes 1 and 2 
in Eqs. (|3-4[) and (|3-5I) . We consider the case that the low-dimensional resonance 
island appears located at Pi = 0, P 2 ^ so that we can easily check the validity of 
our analysis numerically. If one considers the case l|3-lj) . then there is no resonance 
island up to this approximation. Furthermore one considers the cases 1)3- 2 Jl and 
(|3-3|) . there are resonance islands described by 4-dimensional maps. For the case of 
Pi = 0, P 2 ^ 0, a secular term (oc en 2 ) appears in Eq. ()3-4|) . however, does not appear 
in Eq. (|3-5|) . To remove the secular term in Ea. H3-4j) . we define the renormalized 
variable 

a™ := a\ + e—n 2 sin(27rai). (3-6) 

47T 

From the definition (|3-6j) . we obtain the symplectic renormalization map 

La\ = e—^u{2i:a r {). (3-7) 

This procedure to obtain the renormalization map Q3-7|) is the same as the case for 
the standard map in The renormalization map (|3-7|) is symplectic and has the 
form of the standard map. Since the renormalization map of the standard map is 
the standard map itself (see Eas. (|2-l|) and (|2-4|0 . a self-similar structure of the phase 
space could be obtained when one considers the renormalization map Q3-7JI . 

From Eas. (|3-4[) - (|3-5|) and Eq. (|3-6j) . we have the following relation between the 
renormalized variable a" and the original variables (x^x^) 

„ C sin(2ira? + 27r(a2 + nPo)) „, ,, . , , . 

X 1 = a1 + e- ^ 1 \ 2 V > +0(e 2 ), modi 3-8 

47T cos(z7rP 2 j — 1 

( A 2 sin(2^(a 2 + nP 2 )) C sin(2^ + 27r(q 2 + nP 2 )) 
X2 _ a2 + nF2 + e y^ cos{27r p 2) _ 1 + 47r cos(27rP 2 )-l 

+C(e 2 ), (modi) (3-9) 

and the other original variables y™ and y% are obtained as y™ = x™ — x™~ 1 and 
= x 2 ~ x 2~ 1 - The unstable manifold of the renormalization map (|3-7|) and 
Eqs. H3-8|) - ()3-9|) reproduce the low-dimensional resonance island in the ambient 
phase space of the system. Figs. |2] and 01 show that our analysis can predict the 
low-dimensional resonance island even in 4-dimensional phase space. Due to this 
success of the prediction, we understand that the low-dimensional resonance island 
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Fig. 2. The phase portrait near the resonance island characterized by Pi = 0, Pi 7^ 0, parameters 
are P 2 = 0.18(/ Q),eA 1 = 0.04, eA 2 = 0.05, sC = 0.02. The number of the initial condition is 
1. (a): obtained by the renormalization method, and (b): calculated using the original map. 




Fig. 3. The phase portrait near the resonance island characterized by Pi = 0, P2 7^ 0. The 
parameters are given in the caption to Fig. [5] (a): obtained by the renormalization method, 
and (b): calculated using the original map. 



has the same functional form as in the 2-dimensional standard map. To confirm the 
description of the approximate low-dimensional resonant island further, we plot a 
projection of the manifold onto the space spanned by x\ and y\ as shown in Fig. |1J 
and it certainly suggests the success of the analysis with the perturbation strength 
< |e| < 0.01. Although we do not show any figure, our perturbative analysis fails 
when e is bigger than ~ 0.01 for P± = and P2 = 0.18. Taking into account 
Fig. ^ one finds that the regime < |e| < e* in which we can use a perturbation 
analysis for the Froeschle map is narrower than that for the standard map in 33 
A mathematical estimate is needed to determine the actual validity of our analysis 
again, as stated in £j2j On the other hand, according to the expressions for the naive 
perturbation solutions H3-3|H()3-5|) . it is natural to expect that the renormalization 
map (|3-7j) fails when P2 approaches to 0, and we have numerically confirmed that e* 
generally depends on P2 for Pi = 0. A situation where the integration constant P2 
narrows the regime < |e| < £* does not appear in 2-dimensional maps, and is one 
of features in high-dimensional maps. 
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Fig. 4. The phase portrait near the resonance island characterized by Pi = 0, Pi 7^ 0. The 
parameters are given in the caption to Fig. [5] (a): obtained by the renormalization method, 
and (b): calculated using the original map. 



§4. Conclusions 

In this article, we have analyzed a 2- and a 4-dimensional symplectic maps us- 
ing a renormalization method. For both maps, we have successfully reduced their 
resonance island analytically and the perturbative reductions have been confirmed 
numerically. Although a resonance island is the most fundamental object to study 
global Hamiltonian chaos for a system with any degrees of freedom, the island struc- 
ture cannot be visualized for a system with many degrees of freedom. This study 
shows how a low-dimensional island appears in the Froeschle map, at least in a 
nearly integrable regime. The method which we have used here is also useful for a 
wider class of symplectic maps, and then this analysis shed light on the study of the 
Hamiltonian systems with many degrees of freedom. 

Acknowledgments 

The author thanks the members of the laboratory of dynamical systems theory, 
Kyoto University, for their comments and encouragements. The author has been 
supported by a JSPS Fellowship for Young Scientists. 

References 

1) A.J. Lichtenberg and M.A. Lieberman, Regular and Chaotic Dynamics (Springer, Berlin, 
1992). 

2) L.Y. Chen, N. Goldenfeld and Y. Oono, Phys. Rev. Lett. 73 (1994), 1311 ; L.Y. Chen, N. 
Goldenfeld and Y. Oono, Phys. Rev. E. 54 (1996), 376 . 

3) Y. Oono, Int. J. Mod. Phys. B. 14 (2000), 1327 . 

4) T. Kunihiro, Prog. Theor. Phys. 94 (1995), 503 ; T. Kunihiro, Prog. Theor. Phys. 97 
(1997), 179 ; T. Kunihiro and J. Matsukidaira, Phys. Rev. E 57 (1998), 4817. 

5) S. Goto, Y. Masutomi and K. Nozaki, Prog. Theor. Phys. 102 (1999), 471. 

6) Y. Masutomi and K. Nozaki, Physica D. 151 (2001) 44; S. Kawaguchi, Prog. Theor. Phys. 
113 (2005) 687. 

7) S. Goto and K. Nozaki, Prog. Theor. Phys. 105 (2001), 99; S. Goto and K. Nozaki, J. 
Phys. Soc. Jpn. 70 (2001), 49; S. Goto, K. Nozaki and H. Yamada, Prog. Theor. Phys. 
107 (2002), 637; S.I. Tzenov and R.C. Davidson, New Journal of Physics, 5 (2003), 67; 
S. Goto and K. Nozaki, Physica D. 194 (2004), 175; T. Maruo, S. Goto and K. Nozaki, 



8 



Shin-itiro GOTO 



Prog. Theor. Phys. Ill (2004), 463. 

8) D.S. Broomhead and G. Rowlands, J. Phys. A. 16 (1983), 9. 

9) C. Froeschle, Astrophys. SpaceSci. 14 (1971), 110. 

10) For exmaple, B.P. Wood, A.J. Lichtenberg and M.A. Lieberman, Phys. Rev. A 42 (1990) 
5885; S. Honjo and K. Kaneko, Advances in Chemical Physics 130 part B (2005)437. 



